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In this study, slow pyrolysis characteristics of olive oil pomace were investigated by using 
thermogravimetric analysis coupled with mass spectrometry. The major pyrolysis prod¬ 
ucts identified are H 2 , CH 4 , CO, C0 2 . Their evolution profiles with respect to temperature 
enable a realistic evaluation of weight loss results. The chemical structure of the olive oil 
pomace was analyzed using FTIR. Thermogravimetric analysis results have been utilized to 
determine kinetic parameters by using model fitting (Coats Redfern method), model free 
(ASTM E698, Flynn—Wall—Ozawa, Friedman methods) and nonlinear regression analysis 
approaches. Comparative evaluation of the kinetic results indicate that the multivariate 
regression analysis is an appropriate method to derive kinetic models which give reliable 
results for the whole temperature range especially in multi-step reactions. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

Sustainable primary and secondary energy generation from 
biomass is of utmost importance considering the diminishing 
availability of fossil fuels and the higher sensibility towards 
environment preservation from pollutants generated by con¬ 
ventional energy production. Renewable energy derived from 
biomass reduces reliance on fossil fuels and it does not add 
new carbon dioxide to the atmosphere. Biomass is a term for 
all organic material that stems from plants including algae, 
trees and crops that are susceptible to be converted into en¬ 
ergy [1], 

Biomass thermolysis has the potential to reduce our 
dependence on fossil resources by providing alternative fuels 
and specialty chemicals. Biomass pyrolysis is an important 


step in all thermolysis processes such as gasification, lique¬ 
faction and combustion [2], Therefore know-how on pyrolysis 
products and related kinetics is essential in various energy 
and chemical production technologies [3], Pyrolysis is an 
extremely complex process which goes through a series of 
reactions and can be influenced primarily by the temperature 
profile and pressure of the process [4—7]. The differences 
observed in the pyrolysis behavior of the ligno-cellulosic 
samples can be attributed primarily to the variations of 
hemicellulose, lignin and cellulose content of the biomasses 
under investigation [8]. 

Thermogravimetric analysis (TGA) is a commonly applied 
thermo analytical technique in solid-phase thermal degrada¬ 
tion studies [9], and it has gained widespread application in 
thermal studies of biomass pyrolysis [10—15], TGA coupled 
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with mass spectrometry (MS) enables the identification of 
gaseous species formed during pyrolysis. 

Thermogravimetric analysis (TGA) results can be used for 
kinetic analysis of thermolysis processes. In literature, 
numerous works describe TGA based kinetics of different 
materials such as plastics [16—18], inorganics [19—21], urban 
waste [22-24], fossil fuels [25-27] and biomass [4,28,29] during 
thermal degradation. 

A kinetic model is necessary to predict the thermolysis 
behavior, to design the necessary equipment and to deter¬ 
mine their operation conditions [30]. 

This paper presents results related to slow pyrolysis of 
olive oil pomace, the solid residue of olive oil production. This 
residue includes olives (6—8%), water (20—33%), seeds and 
pulps (59—74%) depending on the olive oil extraction pro¬ 
cesses [31]. The annual world production of olive oil is about 3 
million tons (statistics from year 2009) yearly and Turkey has 
ca. 5% share in olive oil production [32], Slow pyrolysis of olive 
oil pomace occurs during torrefaction, active carbon forma¬ 
tion and char gasification. Hence, the results of this study are 
applicable to primary and secondary fuel productions from 
pomace. The produced gas mixture can be converted to syn¬ 
thesis gas which can be catalytically processed to various 
hydrocarbon products. The solid residue can be utilized as 
active carbon for adsorption and catalysis. 


2. Materials and methods 

2.1. Samples 

The olive oil pomace used in this study was supplied by the 
Turkish olive oil producers. 

The dried olive oil pomace samples were grinded with a 
knife mill to get the final size of the sample below 250 p. 
Proximate and ultimate analysis results were carried out 
according to ASTM Standard D5142-04 (TGA 701, LECO). Ul¬ 
timate analysis was conducted using ASTM Standard D5373- 
2 (Truspec CHN-S, LECO). The results are summarized in 
Table 1. 

2.2. Experimental apparatus and procedures 

Thermogravimetric analysis of olive oil pomace was carried 
out using Netzsch STA 409 coupled with AELOS 403C QMS 


Table 1 - Characteristics of the olive oil pomace 
basis). 

(wt%-dry 

Proximate analysis 


Volatile matter 

74.26 

Fixed carbon 

18.95 

Ash 

6.79 

Ultimate analysis 


C 

54.72 

H 

5.24 

O 

31.97 

N 

1.14 


mass spectrometer. A sample mass of 20.0 mg was used in 
each experiment. Argon was used as a carrier gas with a flow 
rate of 45 mL/min for pyrolysis experiments. The olive oil 
pomace samples were heated from 25 °C to 900 °C at four 
selected heating rates, 1.0, 2.5, 5.0 and 10 °C/min. The sample 
mass used is high enough to represent the olive pomace 
sample and low enough to decrease the heat transfer limi¬ 
tation. The applied low heating rates also serve to minimize 
the impact of heat transfer resistance as much as possible 
[33]. 

A portion of the volatiles formed as a result of thermal 
decomposition was transferred to the ion source of the QMS 
403C through a quartz capillary line heated to 250 °C. The ml 
z ratios of selected products were as follows: ml z 2, 15, 28, 
44. 

The samples were mixed with KBr for solid FTIR analysis 
following the standard procedure. Several pomace pellets at 
various mixture ratios were prepared. Results presented here 
refer to a pomace KBr 1:30 mass ratios. 

2.3. Kinetic analysis 

Kinetic analysis is traditionally expected to produce an 
adequate kinetic description of the process in terms of the 
reaction model and the Arrhenius parameters, where fe 0 is the 
pre-exponential factor, f(a) a function called the reaction 
model, E„ the activation energy and R the gas constant. These 
three components ( f(a ), E„, and fe 0 ) are also called the “kinetic 
triplet” [34]. 

The rate of heterogeneous solid-state reactions can be 
described by the Arrhenius equation where t is time, fe(T) is 
the temperature-dependent constant, which describes the 
dependence of the reaction rate on the extent of reaction, a 
(Eqn. (1)). 

J = fe(T)/(«) (1) 

The extent of reaction, during the thermal degradation can 
be defined as the ratio of weight loss at time t to the total 
weight loss corresponding to the degradation process (Eqn. 
(2)). 



J=fe 0 e^/(«) (3) 

There are many methods for analyzing non-isothermal 
solid-state kinetic data from TGA [35,36]. Knowledge of the 
kinetic parameters, such as the reaction rate and activation 
energy, is one of the keys to determine reaction mechanisms 
in solid phase [34]. 

Kinetic analysis techniques can be classified in three 
groups as model-fitting (i.e., identification of a kinetic reaction 
model), iso-conversional (i.e., model-free) and nonlinear 
regression analysis (Fig. 1). 

2.3.1. Model fitting methods 

In the model-fitting method, the f(a) term is determined by 
fitting various reaction models to experimental data. 
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Subsequently, the fe(T) can be evaluated by the chosen form 
of f(a) [34]. The model-fitting methods are widely used 
because of their capacity to directly determine the kinetic 
triplet [34]. In the model fitting technique, the apparent 
activation energy is assumed to remain constant throughout 
the whole temperature range. Equation (3) can be trans¬ 
formed into non-isothermal rate expressions describing re¬ 
action rates as a function of temperature at a constant 
heating rate, (1: 


da 

dT 


>%(«) 


Upon integration 


(4) 


2.3.2.1. ASTM E698. Kissinger’s method is widely used as a 
standard ASTM E698 procedure. Its application in pyrolysis 
may be difficult due to overlapping reactions [48], 

The value of activation energy is calculated from a plot of 
ln(/3/T^) against 1000/T m for different heating rates ((3), where 
T m is the temperature peak of the DTG curve. The equation is 
as follows: 


E„ 

RTm 


(9) 


E a can be calculated from the slope of the plot, which is equal 
to E„/R. 


9 ( “) = t/ B ® dT (5) 

If (E„/RT) is replaced by x and integration limit is 
transformed: 



Equation (6) can be written as: 


9(“) = ^P( x ) ( 7 ) 

p(x) has no analytical solution but has many approximations 
[37,38], one of the most popular is the Coats Redfern method 
[39-42], 

2.3.1.1. Coats Redfern method. This method utilizes the 
asymptotic series expansion for approximating the expo¬ 
nential integral in Eqn. (8), giving 

< s > 

g(a) versus 1/T, gives E a and fe 0 from the slope and intercept, 
respectively. The model that gives the best linear fit is selected 
as the chosen model [43]. 

2.3.2. Model free methods 

Model-free iso-conversional methods allow the activation 
energy to be estimated as a function of a without choosing the 
reaction model. The basic assumption of these methods is 
that the reaction rate for a constant extent of conversion, a, 
depends only on the temperature [44-47]. 



2.3.2.2. Friedman method. Friedman method is a linear dif¬ 
ferential method based on iso-conversional technique which 
assumes the apparent activation energy to remain constant. 
The Friedman method is in general terms as follows: 

*-»(£) = *•'**«> <“> 

It is assumed that the conversion function f(a) remains 
constant, which implies that biomass degradation is inde¬ 
pendent of temperature and depends only on the rate of mass 
loss. A plot ofln[/J(da/dT)] versus 1/T yields a straight line, the 
slope of which corresponds to -E a /R [44], 

2.3.2.3. Flynn Ozaiva methods (FWO). The Flynn-Wall- 
Ozawa (FWO) method [45,46,49,50] is an integral iso- 
conversional technique that assumes the apparent activa¬ 
tion energy to remain constant. 

Apparent activation energy (E„) is obtained from a plot of 
ln(ft), versus 1000/T aji , which represents the linear relation 
with a given value of conversion at different heating rates. 

■"“'-'"(w)- 5 - 331 - 1 - 052 ^: < 12 > 

The subscripts i and a denote the given value of heating rate 
and the given value of conversion, respectively. 

2.3.3. Nonlinear regression analysis 

Non-linear regression (NLR) allows a direct fit of the model to 
the experimental data without a transformation and there are 
no limitations with respect to the complexity of the model. For 
this reason, NLR method can be used and iterative procedures 
can be employed for estimating the reaction model of olive oil 
pomace. 

The calculation of the regression values is carried out by 
means of a fifth-degree Runge—Kutta procedure using the 
Prince-Domand embedding formula for automatic optimiza¬ 
tion of the number of supporting digits [51,52]. To minimize 
the deviance, least squares (LSQ) method is used and a smooth 
convergence is ensured. 

The calculated curves were obtained by means of NLR. 
These curves were fitted to the experimental ones and cor¬ 
rected with LSQ. During this procedure quality kinetic pa¬ 
rameters of initial values are optimized [52]. 


Fig. 1 - Kinetic analysis techniques. 
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Mass/% 



Temperaturerc 


Fig. 2 - Weight versus temperature for various heating 


3. Results and discussion 

3.1. Thermogravimetry 

Olive oil pomace like other biomass sources mainly consists of 
hemicellulose, cellulose and lignin together with some resins 
and minerals. Hemicellulose of olive oil pomace begins to 
decompose after dehydration at lower temperatures (<150 °C) 
followed by depolymerization at higher temperatures [49]. 
Cellulose is decomposed in a way similar to hemicellulose: 
Decomposition dominates at low temperatures (<300 °C) and 
slow heating rates. At higher temperatures, depolymerization 
and fragmentation are dominant. Inorganic species such as 
cations, bases, acids and salts can have a major impact on 
these reactions [53—55], Lignin-the strengthening component 


DTG /(%/min) 
-0 7 • 



100 200 300 400 500 600 


T emperature /°C 

Fig. 4 - Three-component devolatilization model for olive 
oil pomace at a heating rate of 1 °C/min. 


of the cell wall is thermally more stable than cellulose and 
hemicellulose and yields more char and a higher fraction of 
aromatic compounds [53—56]. It is a polymer of hydroxy and 
methoxy substituted propyl phenol units [57], At lower tem¬ 
peratures decomposition dominates, while at higher temper¬ 
atures cracking, reforming, dehydration, condensation, 
polymerization, oxidation and gasification reactions take 
places [58—65]. The pyrolysis of cellulosic components pri¬ 
marily occurs around 350 °C, while that of lignin covers the 
whole temperature range (150-800 °C) the emphasize being in 
the high temperature region [2]. Moderate temperatures of 
around 500 °C and short reaction times of up to 2 s are needed 
to mitigate the secondary vapor phase reactions during 
biomass pyrolysis [4]. However, the effect of secondary vapor 



Fig. 3 - Pyrolysis steps at a heating rate of 1 °C/min. 
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Fig. 5 - Evolution of H 2 (m/z, 2) with temperature increase at a heating rate of 1 °C/min. 


phase reactions on pyrolysis product distribution with ther¬ 
mal analysis is poorly known primarily due to the utilized low 
sample masses. 

TGA results of olive oil pomace pyrolysis under argon at¬ 
mosphere at four different heating rates are given in Fig. 2. As 
the heating rate increases tenfold from 1 °C to 10 °C the TG and 
DTG profiles remain almost the same, so are the gas evalua¬ 
tion profiles. The shift is to slightly higher temperatures with 
increasing heating rates. 

Fig. 3 presents the TG-DTG results for 1 °C/min heating 
rate. For all heating rates, the DTG curve shows two consec¬ 
utive pyrolysis steps; primary pyrolysis of hemicellulosic and 
cellulosic components followed by mainly lignin pyrolysis. 
The main DTG peak corresponds primarily to cellulose 


decomposition; the shoulder at lower temperatures can 
mainly be attributed to hemicellulose devolatilization. Con¬ 
tributions from decomposition of lignin and extractives are 
visible as wide shoulders on the main devolatilization 
domain, due to the very broad decomposition range of these 
components [56]. 

The devolatilization profiles of hemicellulose, cellulose 
and lignin for olive oil pomace can be further explained by 
applying a three-component devolatilization model to the 
DTG peak [65—68]. Fig. 4 presents the model applied to olive oil 
pomace pyrolysis at 1 °C/min. The main DTG peak corre¬ 
sponds primarily to cellulose decomposition; the shoulder at 
lower temperatures can mainly be attributed to hemicellulose 
devolatilization. Contributions from decomposition of lignin 



Fig. 6 — Evolution of CH 4 (m/z, 15) with temperature increase at a heating rate of 1 °C/min. 
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Fig. 7 - Evolution of CO (mlz, 28) with temperature increase at a heating rate of 1 °C/min. 


and extractives are visible as wide shoulders on the main 
devolatilization domain [65—68], 

As Fig. 4 indicates hemicellulose of olive oil pomace is 
easily degraded, and its pyrolysis occurs around 150-295 °C. 
The pyrolysis of cellulose mainly happens around 245—335 °C, 
while that of lignin covers the whole temperature range 
(150-480 °C). 

Hence, the first lower temperature step is associated with 
the decomposition of hemicellulosic and cellulosic materials; 
this step also contains part of the low temperature lignin 
decomposition products. The global maximum degradation 
rates are indicated by the DTG peaks. The second higher 
temperature decomposition step is less pronounced then the 


lower temperature step and indicates primarily lignin 
pyrolysis. 

During the primary pyrolysis of olive oil pomace thermal 
cracking of chemical bonds in the individual constituents of 
olive oil pomace, which are mainly cellulose, hemicelluloses 
and part of the lignin, takes place and permanent gas species 
(e.g. C0 2 , CO, CH 4 ) and condensable species at ambient con¬ 
ditions (several organic compounds and water) are evolved. 

The secondary pyrolysis of olive oil pomace takes place in 
the range of 335—480 °C along with the formation of a solid 
char residue. During the secondary pyrolysis, serial and par¬ 
allel reactions take place, occurring either heterogeneously 
or homogeneously, as for example cracking, reforming, 


TG/% 


DTG /(%/min) 
Ion Current *10-« /A 



Fig. 8 - Evolution of C0 2 (m/z, 44) with temperature increase at a heating 


of 1 °C/min. 
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Wavenumber cm' 1 


Fig. 9 - The FTIR spectrum of the olive oil pomace. 


dehydration, condensation, polymerization, oxidation and 
gasification reactions [58—65], 

3.2. Mass spectrometry 

Combining TGA with gas-analytical techniques enhances the 
possibilities for correctly interpreting the mechanism of 
thermally induced reactions, which involve the formation of 
gaseous species. It offers the ability to solve a wide range of 
problems and applications based on the determination of the 
composition of complex gas mixtures produced from the 
decomposition of different types of sample. 

Numerical results presented here refer to 1 °C/min heating 
rate. The peak temperatures of the evolved gases increase 
with increasing heating rates. The maximum difference has 
been observed with hydrogen. A tenfold increase of the 
heating rate from 1 to 10 °C/min increased the peak temper¬ 
ature 15%. 

The evolution of H 2 shows a first major peak at 305 °C and a 
shoulder referring to a second peak around 400 °C (Fig. 5). The 
release of H 2 at higher temperature is due to the cracking of 
heavy hydrocarbons (secondary pyrolysis), while the release 
at the lower temperature is due to primary pyrolysis [69], 
Hence, both primary and secondary pyrolysis reactions pro¬ 
duce hydrogen. 

CH 4 evolution occurs throughout the whole pyrolysis 
temperature range (Fig. 6). In the lower temperature region the 
shoulder indicates hemicellulosic CH 4 followed by the peeking 
CH 4 evolution primarily from cellulosic components. CH 4 is 
mainly produced by the cracking of methoxyl groups of the 
lignin at higher temperatures [56,69], 

The evolution of CO commences at around 245 °C, reaches 
its maximum at 335 °C, and then it is prolonged (Fig. 7). The 
release of CO and C0 2 between 245 and 335 °C is probably due 
to the breaking of C—C and C—O bonds which essentially 
originate from hemi-cellulose and cellulose and the release of 


CO at higher temperatures is due to lignin, possibly through 
the release of COOH groups and the breaking of C—O groups 
[62]. Cellulose gives the highest CO yield, due to its high 
carbonyl content [56,69,70], 

C0 2 release is mainly due to the primary pyrolysis, 
while secondary pyrolysis is the main source for the 
release of CO and CH 4 (Fig. 8). Hemicellulose displays the 
highest C0 2 yield because of the higher carboxyl content 
[56,69,70]. 

3.3. FTIR analysis 

The FTIR spectrum of the olive oil pomace is shown in Fig. 9 
olive oil pomace contains primarily alkenes, esters, aro¬ 
matics, ketones and alcohols with different oxygen- 
containing functional groups [56], The FTIR signals with the 
assigned possible compounds and the typical functional 
groups are listed in Table 2 for a reference [2,71,72], 


1 Table2-The 


i functional groups 

of olive oil pomace. 1 

Wave number 


Functional 

Compounds 

(cm- 1 ) 


groups 


3273 


O-H stretch 

Acid; ether linkage 

2927 


C-H stretch 

Alkyl 

2360 



Linked water 

1649 


CH 2 bend 

Alkanes 

1515 


O-H bend 

Acid 

1378 


C-H bend 

Alkyl 

1321 



Xylan 

1249 


C—O—C stretch 

Pyranose ring skeletal 



vibration 


1039 


C-0 stretch 

Hydroxyl group 

665 


C-H bend 

P-D-Xylose 










79 


175 



3.4. Kinetic results 
3.4.1. Model fitting method 

The activation energy (E a ) and the pre-exponential factor (fe 0 ) 
were obtained assuming the olive oil pomace decomposition 
as a first order reaction with the Coats Redfern method uti¬ 
lizing the TG data at 1 °C/min. 

The first temperature range of interest is chosen as 
250 °C—310 °C based on the DTG results. The low temperature 
limit corresponds to the maximum decomposition rate of 
hemicellulose. The temperature range includes the maximum 
decomposition rate of cellulose at 292 °C and ends when cel¬ 
lulose decomposition is completed (Fig. 4). Fig. 10 shows a plot 
of ln( g(a))-21n(l/T) versus 1/T for pyrolysis kinetics of olive oil 
pomace between 250 °C—310 °C. The results indicate an acti¬ 
vation energy of 179.15 kj/mol and a pre-exponential factor of 
3.09 x 1012 s h 

In the second olive oil pomace pyrolysis kinetics calcula¬ 
tion approach, two consecutive temperature ranges which 



pyrolysis kinetics of olive oil pomace between 
150 °C—335 °C. 



Fig. 12 - Phase 2: A plot of ln(g(«))-21n(l/T) versus 1/T for 
pyrolysis kinetics of olive oil pomace between 
335 °C—480 °C. 


correspond to primary and secondary pyrolysis phases were 
selected. Figs. 11 and 12 show the associated kinetic calcula¬ 
tion results. The lower temperature range covers 
150 °C-335 °C (phase I) while the higher temperature range is 
between 335 °C—480 °C (phase II) .The values of the activation 
energy and pre-exponential factors were calculated as 
78.70 kj/mol and 79.93 kj/mol and a pre-exponential factor of 
7.60 x 102 s 1 and 2.87 x 10 s 1 for pyrolysis phases I and II, 
respectively. 

3.4.2. Model free methods 

A series of non-isothermal TG measurements with different 
heating rates result in a set of data which enables the collec¬ 
tion of data for the same degree of conversion (a) at different 
temperatures [34]. Based on this idea, the activation energy 
(E a ) and pre-exponential factor (fe 0 ) are obtained using Kis¬ 
singer, Friedman and Flynn—Wall—Ozawa (FWO) methods. 
The data (/?, a, T) for kinetic calculations explained in Section 
2.3 are taken from the TG curves. 


log (Healing rate/(K/min)) 



Fig. 13 - A plot of log (/3) against 1000 K/T m at four different 
heating rates (1, 2.5, 5,10 °C/min). 
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Fig. 14 - The plot of log (da/dt) against 1000 K/T at four 
different heating rates (1, 2.5, 5, 10 °C). 


3.4.2.1. ASTM 698. Kinetic parameters are obtained from TG 
curves at four different heating rates (1, 2.5, 5,10 °C per min¬ 
ute) using ASTM E698. The slopes and the intercepts of the 
straight lines (Fig. 13) give the activation energy and the pre¬ 
exponential factor as 170.99 kj mol -1 and 6.46 x 10 12 s \ 
respectively. 

3.4.2.2. Friedman analysis. The graphs ln(da/dt) versus 1000/T 
yield straight lines for each conversion a with slopes directly 
proportional to the reaction activation energy (E a ) (Fig. 14). 

The activation energy and pre-exponential factors remain 
almost constant around 173 kj/mol and 1.65 x 10 14 s 1 in the 
conversion range of 0.05 < a < 0.65. An increase of these 
values is observed at higher conversions. Fig. 15 shows these 
lines at different a and the calculated results using both Eqs. 
(10) and (11). The upper and lower curves in Fig. 15 indicate the 
statistical errors associated with Friedman analysis for pre¬ 
exponential factors and activation energies, respectively, at 





Fig. 16 - The plots of log(/3) against 1000 K/T at four 
different heating rates (1, 2.5, 5,10 °C). 


fractional mass losses. Our TG results suggest that results are 
acceptable up to 0.7 conversions. Beyond this range (a > 0.70), 
these values start to fluctuate. 

3.4.2.3. Flynn—Wall—Ozawa. The graphs ln(/3) versus 1/T 
and ln(da/dt) versus 1/T both show straight lines (Fig. 16) with 
slopes directly proportional to the reaction activation 
energy (E„). The values of the activation energy and pre¬ 
exponential factors are increasing slowly in conversion 
range of 0.05 < a < 0.65, and beyond this range (a > 0.65), both 
kinetic values increase rapidly. The upper and lower curves in 
Fig. 17 indicate the statistical errors associated with 
Flynn—Wall—Ozawa analysis for pre-exponential factors and 
activation energies, respectively, at fractional mass losses. 

3.4.3. Non-linear regression analysis 

Multivariate non-linear regression (NLR) program was applied 
to determine the mechanism and the corresponding kinetic 
parameters of slow pyrolysis. Simulations were performed by 
using “Netzsch Thermokinetics: A Software Module for the 
Kinetic analysis of Thermal Measurements”. The non- 


E/(kJ/mol) log(WsM) 





Fig. 15 — Friedman analysis-energy plot, /3 “1, 2.5, 5 and 
10 °C/min”. 


E/(kJ/mol) k>g(k./sM) 



FractMass Loss 


Fig. 17 - FWO analysis-energy plot, /3 “1, 2.5, 5 and 10 °C/ 
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isothermal data recorded at four heating rates were brought 
together during analysis and the relevant differential equa¬ 
tions of the reaction rates were numerically solved and the 
kinetic parameters were iteratively optimized. The calculated 
curves were obtained by means of non-linear regression. 
These curves were fitted to the experimental ones and cor¬ 
rected with least square method. During this procedure initial 
values are optimized in order to get high fitting quality kinetic 
parameters. The calculations were performed for conversions 
of 0.05 < a < 0.70. The proposed four steps reaction model is 
given in Table 3. For all cases the reaction order is changed 
between 0.5—2.5 to find the best fit. 

Fig. 18 shows the fit of the values calculated with non¬ 
linear regression analysis to the experimental TG data at the 
four selected heating rates. The correlation coefficient is 
0.9977. 

Table 4 summarizes the activation energy and pre¬ 
exponential factor results for the before mentioned kinetic 
methods along with the non-linear regression method. 


4. Concluding remarks 

Investigation of the slow pyrolysis kinetics of olive oil pomace 
using thermo-gravimetric analysis coupled with mass spec¬ 
trometry enables co-current determination of the pyrolysis 
based weight loss and the evolution of selected gaseous 



products as a function of temperature. Both weight loss and 
gas evolution (in this case H 2 , CF1 4 , CO, C0 2 ) profiles are useful 
in the interpretation of the kinetic results. In other words, 
TGA/MS experiments can be designed to understand the ki¬ 
netics of thermolysis processes. 

Pyrolysis reactions are not simple one-step processes; they 
constitute a combination of serial and parallel steps. Hence, 
kinetic parameters change during the course of the reactions. 
The most applied kinetic method in literature is the Coats 
Redfem method. Our studies indicate that the temperature 
range of interest is of utmost importance for the results of this 
method. The activation energy and pre-exponential factor 
results with Coats Redfern based on maximum conversion are 
in good agreement with the model free methods. Kinetic re¬ 
sults obtained with Coats Redfem for two wide pyrolysis 
temperature ranges are much lower. The results of the chosen 
four steps kinetic model are also in the range of the former 
results. Only model free methods estimate the potential 
model(s) according to which the reactions occur. All ap¬ 
proaches give information about the pyrolysis kinetics. One 
should first of all comparatively evaluate them and then the 
method most suitable for the research purpose should be 
selected. 


Mass/% 



Temperature/°C 


Fig. 18 - The calculated values to the experimental data at the four selected heating rates. 
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Nomenclature 


a the given value of conversion 

E a apparent activation energy (kj mol 1 2 ) 

f(a) reaction model (function expressing the dependence 
of the reaction rate on the conversion) 
g(a ) integrated reaction model 

i denote the given value of heating rate 

fe 0 pre-exponential factor (s 

fe(T) temperature-dependent rate constant {s 1 ) 

p(x) temperature integral 

R universal gas constant (8.3144 x 10 3 4 5 6 7 8 9 10 kj mol 1 K 1 ) 

t time (s) 

T temperature (K) 

x E„/RT 

Greek letters 

a conversion 

/3 heating rate (°C/s) 
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